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Abstract. We show that the semiclassical approach to chaotic quantum 
transport in the presence of time-reversal symmetry can be described by a 
matrix model, i.e. a matrix integral whose perturbative expansion satisfies 
the semiclassical diagrammatic rules for the calculation of transport statistics. 
This approach leads very naturally to the semiclassical derivation of universal 
predictions from random matrix theory. 
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1. Introduction 

We consider wave scattering through a system in which the corresponding ray 
dynamics is strongly chaotic. Under a minimal information statistical approach, it is 
natural to model the scattering S matrix as a random matrix, which must be unitary 
by conservation principles [lElE]- The quantity Sij is the scattering amplitude from 
channel j to channel i. In the presence of time-reversal symmetry, scattering from i 
to j is equivalent to scattering from j to i and so S must also be symmetric. The 
ensemble of unitary symmetric complex matrices has a natural probability measure 
on it, and is known and the Circular Orthogonal Ensemble (COE) of random matrix 
theory (RMT). 

RMT predicts successfully several scattering observables, in agreement with 
numerical simulations and experimental results, such as average conductance, 
conductance fluctuations, average shot-noise and higher counting statistics. This 
approach has been reviewed in jl], and some recent results include MB [30 ED]. 
It is rather flexible, and can be adapted in order to treat the statistics of time delay 
[HE3E3!, and to consider the presence of superconductors [UlCESj, non-ideal contacts 
EZl E3 El EDI El], graphene [Ml |23], etc. 

Recovering RMT results from chaotic scattering trajectories has long been a 
central problem for the semiclassical approach to quantum mechanics, in which Sij 
is expressed as a sum over paths leading from j to i [Ml [2S]. In order to reproduce 
quantum effects, it is necessary to consider trajectories that are action-correlated on 
the scale of h. This started to be done perturbatively in i26| [27] and was shown to give 
the exact result for the simplest observables in [H J2D1 |3D]. Attention then turned to 
more general transport statistics [SI El 133 El, until complete equivalence was shown 
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between semiclassics and RMT [351 ' 5HT] (in the meantime, semiclassics was able to go 
beyond RMT, incorporating effects due to finite Ehrenfest time, see e.g. E3 ESI ESI)- 

These semiclassical works rely on elaborate manipulations with diagrams and/or 
permutations. A more direct demonstration of the RMT-semiclassics equivalence 
was presented in [40], based on a matrix model formulation of the semiclassical 
approximation. This is a matrix integral whose diagrammatic expansion satisfies 
exactly the same diagrammatic rules as the semiclassical calculation of transport 
observables, and which turns out to be equivalent to usual RMT. 

However, the treatment in m concerns only systems where time-reversal 
symmetry is broken. The purpose of this work is to extend the matrix model approach 
to time-reversal symmetric systems, establishing the RMT-semiclassics equivalence in 
a direct way for this universality class. 


2. Usual RMT treatment 


Suppose a chaotic cavity, coupled to two ideal leads supporting N± and open 
channels, having a M— dimensional S matrix where M = (N\ + Nf). The energy 
of the incoming wave is E, and the classical dynamics in the cavity is assumed fully 
chaotic at this energy. Under this assumption, RMT assumes the statistical properties 
of transport to be independent of E. 

We shall consider the quantities 

-P(b j) = <5iii2«5i 3 i4 • • • Sii n - 1 ii n Sj 1 j 2 Sj a j i ■ ■ ■ S j2n _ ij2n , (1) 

which can be used to expand any generic observable (all matrices are taken at the 
same energy). For instance, the (dimensionless) conductance and shot-noise are given 

by 


9 =Y S ‘ 


c* 

22 5 


p = 9- Y 


c. . c* c. . c* 

<-^122 ^ 3 * 2 ^ 3 ^ 4 ^ 21^4 5 


( 2 ) 


where the sums over Z2,«4 (*i, *3) run over the N\ incoming (IV2 outgoing) channels. 

According to the random matrix theory approach, for time-reversal invariant 
chaotic systems S is uniformly distributed in the COE(M). The interest then lies in 
the average value of m over this ensemble. The simplest case is 


(I>%| 2 )coe(m) - M Y 


(3) 


which leads to the well known prediction N 1 N 2 KM + 1) for the average conductance. 
In general, the average vanishes unless j is equal to some permutation of i. For 
instance, 


(*^12 ^12 ^34*^34) COE (M) 


M + 2 

M(M + 1)(M + 3) ’ 


(4) 


and 


(<5>12«S 14^34^23 )cOE(M) 


— 1 

M (M + 1) (M + 3)' 


(5) 


There may be several possibilities for the permutation relating j and i. and the 
general expression for our average contains a sum, 



COE(M) 


7r£©2 n 


( 6 ) 
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where & n is the permutation group of n symbols and 

2 n 

= n ( 7 ) 

fc=l 

The quantity Wm(tt), called the Weingarten function of the COE, can be computed 
from explicit character-theoretic formulas [41] or via recurrence relations [42] . 

Since the 5-matrix is symmetric, the value of ([6]) must be invariant under the 
action of any permutation that interchanges j 2 fc-i with j- 2 k ; or any permutation that 
interchanges simultaneously j 2 k-i with j 2 r -i and j 2 fc with j 2 r - The set of all such 
permutations is called the hyperoctaliedral group S) n C & 2 n- For n = 2 it consists of 
the permutations {1, (12), (34), (12)(34)}. In general, fj n has 2 n n\ elements. Thus, if 
£ £ f) n then 

Wm«)=Wm(7t). (8) 

3. Semiclassical Diagrammatics 

The semiclassical approximation to Sij requires trajectories starting at channel j and 
ending at channel i. In the semiclassical expression for m we end up with some 
trajectories ( direct ones) going from * 2 * to * 2 fc-i and some other trajectories ( partner 
ones) going from j 2 fc to j 2 fe-i- Consider now the average value of m over a certain 
energy window, ( P(i,j))E , this window being small in the classical scale but large in 
the quantum scale; as ft —»• 0, constructive interference is required and the result is 
determined by correlations: partner trajectories must have almost the same collective 
action as direct ones. 

The theory of correlated chaotic trajectories has been discussed in detail in a 
number of papers [251 [291 ES3 [321 [36]. Trajectories from correlated sets may differ 
only in small regions (called encounters) in which the direct ones run nearly parallel 
or anti-parallel, while the partner ones have crossings. This ensures they have almost 
the same collective action. In particular, this implies that i and j must be equal up 
to a permutation, a condition already met in the RMT treatment. 

These trajectory multiplets are usually represented by diagrams, in which the 
complicated pieces of chaotic trajectories in-between encounters are depicted as 
simple links. Calculation of any given ( P(i,j))E requires constructing all possible 
contributing diagrams. Most importantly, a diagrammatic rule has been devised for 
the value of any diagram: it is (— l) v M V ~ L , where V is the number of encounters 
and L is the number of links. 

We show two examples in Figure 1. Panel a) shows a contribution to 
5,;-i ,> 2 5* i2 S VJ u iS* u . containing two 2-encounters. One of the encounters involves 
two different trajectories, while the other one is of a trajectory with itself. Notice 
how the direct and partner trajectories run in opposite senses in one of the 
regions. This contribution is only possible because we are assuming the system 
has time-reversal symmetry. Its value is 1/M 4 . Panel b) shows a contribution to 
^ 1 * 2 ^* 312 <^* 5 * 6 ^* 146 > containing a single 3-encounter. This diagram does not 
require time-reversal symmetry, and its value is — 1/M 5 . 

A different diagrammatic representation, more convenient, of correlated 
trajectories uses ribbon graphs l36j . In this case we turn every encounter into a 
vertex, and trajectories are depicted as edges of ribbons. It is sometimes necessary to 
perform twists on some of the ribbons. The diagrams from Figure la,b are represented 
as ribbon graphs in Figure 2a,b. 
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Figure 1. Two examples of semiclassical diagrams required in the calculation of 
transport observables. Direct trajectories (solid lines) are correlated with partner 
trajectories (dashed lines), by means of 2-encounters (left) and a 3-encounter 
(right). This is very simplified: actual trajectories are long and chaotic. 




Figure 2. A different representation of the semiclassical diagrams of Figure 1, 
in terms of ribbon graphs. Encounters become vertices and trajectories become 
edges of ribbons. Notice how one of the ribbons must be twisted in the first case. 


4. Matrix model for the semiclassical approach 

4.1. Wick’s rule 
Suppose the matrix integral 

(f(S))s = \j dSe~^ ssT f(S), (9) 

where M is a parameter and 

Z = jdSe~^ ssT ( 10 ) 

is a normalization factor. This integral runs over general TV-dimensional real matrices, 
without constraints (this is known as the real Ginibre ensemble [43]). Notice that the 
dimension N is not related to channel numbers; we continue with M = N\ + N%. 
Gaussian integrals can be performed exactly, leading to 

(SM s = s -^i. ( 11 ) 

When there are many matrix elements being integrated (the number must be even, 
otherwise the result vanishes), we may use the well known Wick’s rule, 

®<r(2fc —1) G(2fc-1) ^O. a (2k)b a {2k) )s 


2 n \ n 

ns-.*) = e n 

fc=i / s ere®fc=i 


( 12 ) 
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Figure 3. Diagrammatical representation of the average value of 
Tr(SS T ) 2 S ai b 1 S a2 b 2 , and the possible connections between lines, according to 
Wick’s rule. In the upper line a direct connection is made between S ai b 1 and 
S a2 b 2 • There are in total 15 connections, and we show how many times each 
topology appears, due to trivial symmetries. 


where the sum is over all possible matchings among 2 n elements. For example, at 
n = 2 there are three possible matchings, 


cl i 6j_ ^'cl 262 ^CL3 63 ^0464 ^ {*^016^ * 5 'cl 2 b<2 ^ ^^0363 *^0464 ^ 

+ (■^ai&i$a3 63) (*^02 62 ^04 £>4) “I - (^ai61 $0464) (*$0262 ^2363)- 


(13) 


Elements of Wl n can be represented by permutations acting on the trivial matching 
{{1, 2}, {3,4},...}. In this sense, the above matchings correspond to the identity 
permutation, to the transposition (23) and to the cycle (243). A permutation a 
represents a matching if and only if it satisfies er(2z — 1) < cr(2z) and a( 2 i — 1) < 
£7(2 i + 1). Any element of 62 n can be decomposed uniquely as the product of a 
member of Wl n and a member of Sj„, i.e. the set 9Jl n can be seen as the coset & 2 n/?)n- 


4-2. Diagrammatics 


Introduce the following diagrammatical representation to the calculation of (fT^l) . Each 
matrix element S a b is represented by a ribbon, having one edge associated with a 
and depicted with a solid line, and the other edge associated with b and depicted 
with a dashed line. Wick’s rule then tells us to draw all possible connections among 
these lines, and associate to each connection a factor 1/M. When two ribbons are 
connected, their solid edges merge and must have the same label, and likewise for the 
dashed edges. 

For example, consider the average value 


(Tr(SS T ) 2 S aihl S Q2b2 ) 


N 

E 

03046364=1 


(*^0161 *^2262*^0363^0463*^0464*^2364) • 


(14) 


The presence of the trace requires repeated indices among the matrix elements. This 
is taken into account by arranging them around a vertex. We show in Figure 3 the 
ribbons associated with m, and the possible connections in the application of Wick’s 
rule. There are six distinct topologies, each appearing with some multiplicity. The 
calculation results in the value 

(21V 3 + V 2 + 81V + 4)(5 aia2 (i blb2 . (15) 

The power of N counts the number of indices being freely summed over, associated 
with closed lines in the diagrams. 
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This kind of diagrammatics has been extensively discussed in the literature. The 
most common situation is when complex hermitian matrices are used |44| . in which 
case no twists are required in the ribbons (they are orientable) and no distinction is 
made between the edges. Complex non-hermitian matrices 1451 maintain orientability 
but allow different edges. For real symmetric matrices the ribbons are not necessarily 
orientable @61147] and edges are indistinguishable. Our case (general real matrices) 
does not require orientability, and has distinguished edges. 


4-3. Introducing channels 


We are using the labels a and b to identify the ‘trajectories’ of our matrix model. 
We must still introduce the channel labels. To that end, let Q = 1 m © On-m be an 
orthogonal projector, having the iff—dimensional identity in its upper-left corner and 
all other entries equal to zero. Let us also define the matrix 

R = XQSQX\ (16) 


where X is some complex M x M matrix. This leads to 

2 n 2 n M 


II 

k =1 


U Xi k a k S ak b k Xi kbk 

k= 1 a*. ,bfc=1 


(17) 


The indices ik represent the channels. The idea is that, upon taking the average 
over the S matrices, the a-labels will produce the direct trajectories connecting the 
channels, while the 6-labels will produce the partner trajectories, which may connect 
the channels in a different way but will necessarily be correlated with the direct ones. 

For future reference, let us also define the matrix Z = XX T . Notice that this is 
a complex symmetric matrix. The elements of Z will be related to direct trajectories, 
and the elements of Z* will be related to the partner trajectories. 

Before going into further details of this method, let us show it in action for the 
simplest example. 


4 . 4 . The conductance 


Conductance requires the average value of 5,;, , 2 5* i2 . We start by writing 


(R 


2121 -^22 22 


M 

'>s = ^2 

a 1 a 2 b 1 b 2 =l 


Xi 1 aiXi 2a2 Xi i b 1 X,i 2 b 2 (S ai ] Dl S a2 b 2 }g 


) S M Ziii 2 1*2 • 


(18) 


The coefficient 1/M is related to the leading order approximation to Eq.(|3]). 

Next, we can produce diagrams with encounters by including traces of powers of 
SS T . The simplest such case is 

'r(SS T ) 2 R ilil R i2i2 ^ = -±(4+8N+N 2 +2N 3 )Z ili2 Z* ii2 .(19) 

The diagrammatics of this average is the same as the one shown in Figure 3. We have 
multiplied by —M because we already know that each encounter must be accompanied 
by such factor in the semiclassical diagrammatics. We have divided by 4 because of 
the rotation symmetry of the 2-encounter. 

As we have seen, closed lines produce powers of N (the S'-matrices inside the 
trace are not truncated). In a semiclassical interpretation, such closed lines would 
represent periodic orbits, forever trapped inside the system. The true semiclassical 
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diagrams contain only scattering trajectories, and no periodic orbits. We can get rid 
of such orbits by the trick of letting TV — > 0. If we perform this trick, we get Z ixia Zf i2 
with a coefficient of — 1/M 2 , which is the second order approximation to Eq.([3]). 

A triple encounter corresponds to 

( s S T ) 3R nil R l2l2 ^ s = j±Z ili2 Z* ii2 , (20) 

where we have discounted a 6-fold rotation symmetry. The situation with two single 
encounters correspond to 


lim 

TV —>0 


1 

2 


—M 


Tt(SS t ) 2 


Rilil 


_ 7 . . 7 * 

M 3 


( 21 ) 


\ L J /s 

Here there is an extra denominator of 2 to account for the exchange symmetry between 
the encounters. The values of (l20l) and (12T1) together produce a coefficient to Z^Z* t 


which is 1/M 3 , related precisely to the third order approximation to Eq.®. 

In order to produce all possible encounters, we must include all possible traces, 
each one of them multiplied by —M, as required by the semiclassical diagrammatic 
rules. We must also discount an overcounting of 2 q from rotation symmetry around 
vertices of valence q , and take into account the exchange symmetry between encounters 
of the same valence. Fortunately, all this is automatically implemented by means of 
an exponential function. 

Therefore, in analogy to |40j . our semiclassical matrix model for the average value 
of iSq i 2 S* li2 is to compute 


lim ( e 
JV-X) 


-M£„> 2 2yTr[(SS T n 


Rilil R12I2 


( 22 ) 


and to extract the coefficient of Zi^Z* t . In order to compute the conductance, we 
just sum over 


f.5. Higher Moments 


To treat conductance, which involves two channels, we considered the quantity 
Ri 1 i 1 Ri 2 i 2 ■ In order to be able to obtain a semiclassical matrix model for higher 
transport moments, involving n channels, we must consider more general quantities 

of the form Ilfcii Ri k i k - 

As an example, let us look at the diagrammatics of this average value for n = 2. 
Using Wick’s rule, it is easy to see that this is given by 


(-Rilil R%2l2 R^3i3 R>4 74 

I y 7 * 7 y* 

' ^ 1*3 ^ 2123 ^ 2*4 ^2224 


) S ~ M 2 17 * 1*2 ■^ 112 ^* 3 * 4^1314 


7 7 * 7 7 * 1 

^2124^124^2223^2223] 


(23) 


The channel labels appear in different arrangements on the right hand side, 
corresponding to the semiclassical diagrams in Figure 4a. This quantity is thus 
producing all the leading order transport diagrams with four channels simultaneously. 

If we wish to know the leading order approximation to the average value of 
a certain quantity, say S,^ ,; 2 Sf ? - 2 S 1:i ,; 4 S* 3 j 4 . we just extract from the above result 
the coefficient of [Z. ni , 2 Z* i j2 Z i:i i 4 Z* ij4 ] (when performing this coefficient extraction 
procedure, we consider the indices to be independent variables, i.e. we do not care 
for possible coincidences among the elements of i). We shall use the notation [x\f to 
denote the coefficient of x in /. 
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a) 



Figure 4. a) Diagrams related to the leading order approximation to 
(-frii*! Ri 2 i 2 -^* 3 ^ 3 ^ 4 * 4 ) 5 - b) The first two diagrams are responsible for Ea. (|28() . 
while the last diagram is responsible for Eq. m- 


For example, looking at Ea. (l23l) it is clear that 

[Zi^Z i3ii ] / f[ R ikik \ = E Z* ii2 Z* 3ii . (24) 

\fc=l / S 

On the other hand, when we consider \Z* . Z* • 1 Z* , Z,* ,• , it follows that 7 and i 
can only differ by a liyperoctahedral permutation: either we have the set identities 
{ji,j2} = {ii,h}, {j3, J4} = {*3, *4} or the set identities {ji, J 2 } = {h,h},{h, 3 i} = 
{* 1 , 12 }. That is, we have that 

[ z hAt\ K»Ku = E Mm). (25) 


In conclusion, our semiclassical prediction for the leading order approximation to the 
average value of S ik i 2 S* lj 2 S i 3 ii SJ 3j is M ~ 2 Mb/)• 

More concretely, we have, for instance, 

[Z 12 Z* 12 Z 34 Z; A ] (RnR 2 2R33R44)s = fp, (26) 

which is indeed the leading order approximation to Eq.©, and 

[Z 12 Z* 14 Z 34 ZZ 3 ] {R 11 R 22 R 33 R 44 ) s = 0, (27) 

reflecting the fact that there are no semiclassical diagrams contributing at leading 
order to the calculation of Eq.©, and that indeed Eq.© is of order M~ 3 . 

We can produce perturbative corrections, i.e. diagrams with encounters, by 
including traces, as we did for the conductance. For instance, a single encounter 
is introduced as 

lim[Z 12 Z-r 2 Z34Z* 4 ] / EEri (SS T f R kk \ = E (28) 

\ fe=i / s 

which is the first correction to Eq. © , and 

hm [Z 12 Z* 1A Z 34 Z* 23 \ ( ^Tri (SS T f R kk ) = E (29) 

\ fc=1 / s 

which is the leading order contribution to Eq.©. The diagrams corresponding to the 
above results can be seen in Figure 4b. As already discussed, the limit N —» 0 is 
necessary to discard diagrams having periodic orbits. 
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As we saw previously, in order to produce all possible encounters we must 
introduce the factor 2 ^ Tl I(' s ' s For convenience, let us define 


~ IJ Z l2k-1 


k—1 


Then, our semiclassical matrix model is 


P(i, j) ) = lim 

v ' E N-tO L 


3 


G(i), 


where 


2 n 


G(i)= 


k=l / S 


(30) 


(31) 


(32) 


This matrix integral has, by construction, exactly the same diagrammatic formulation 
as the semiclassical approach. 


5. Equivalence with RMT 

We wish to show equivalence between semiclassical theory and random matrix theory, 
i.e. between Eqs. and Eq. ©■ 

We start by noticing that, because of the identity 


det(l-A) = e Trlog(1 - j4) = e - E *>^ Tr(j4 ' ,) , 


we can write 


G(i) = 


dS 


2 n 


det(l - SS T ) M/2 Ri, 


(33) 


(34) 


fc =i 


Next, from the theory of truncated orthogonal matrices it can be established (see 
Appendix) that 

/ 7 o / 2 n \ 

— det(l - SS T ) M G n Sa kbk = ( I] °^ k ) , (35) 

k =1 \k=l / 0(M+ 1) 

provided 1 < a^, bk < M and where the average on the right-hand-side is over random 
matrices uniformly distributed (with respect to Haar measure) in the orthogonal group 
of (M + l)-dimensional real matrices satisfying OO t = 1. 

It is known [3SJ 05] that this orthogonal group average is given by a double sum 
over the set of matchings, 

fi^ = A CT (a)A T (6)yV M +i(o-” 1 r), (36) 

\fc=l / 0(M+ 1) cr,rGOT„ 

where the quantity Wm+i is the Weingarten function of 0[M + 1) and 


A cr(a) = 5, 


a cr(2fc-l) ? a cr(2fc) * 


(37) 


fe=l 


Noticing that 


M 


A <r( a ) Xi k a k — J^[(AA )i a( 2 k _ 1 ) ,i a( 2 k) ~ %c 


'(*) 


ai ,...a2n=l 


k= 1 


k =1 


(38) 
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leads to 

Um G(t)= Wm+i 

a,T£DJl n 


(39) 


We now extract the coefficient of 


Z r Zt 


. The only matching er which survives 


this operation is the identity. On the other hand, we have already seen that 


Z% 


z *(i) = 12 




Thus, we get 


lim 

N —>0 L 


Z 7 Z% 


G(j) = Wm+i(t)S p (t(i)J). 


(40) 


(41) 


We may use the invariance of the Weingarten function VV under the action of the 
hyperoctaliedral, Wm+i(pt) = Wm+ i(t), to group both sums into a single sum over 
the whole permutation group. Defining n = pr, we have 

( p (m)) £ = Wm+i(t r)6 n (ij). (42) 

7T 


It was shown in m that the Weingarten functions of the orthogonal group and 
of the COE are related by the simple identity 


Wm+i(tt) = Wm(tt)- 


(43) 


We are hereby showing that this equality, rooted in the fact that COE(M) can be seen 
as the quotient space U(M)/0(M ), where U(M) is the unitary group, is in fact the 
key to the equivalence between the semiclassical and RMT approaches to quantum 
chaotic transport in the presence of time-reversal symmetry, because it leads precisely 
to the fact that 


( p &j)) E = ( P (U)} 


COE(M) ' 


(44) 


6. Conclusions 

We have extended the matrix model approach to semiclassical quantum chaotic 
transport in order to treat systems with time-reversal symmetry. The hardest part of 
the method is designing a matrix integral with the correct semiclassical diagrammatic 
rules. Once this is in place, it leads in a very direct way to the equivalence to random 
matrix theory. 

This approach may open the way to semiclassical calculations that were previously 
unavailable, and may even provide results beyond RMT. For example, it may be 
adapted to treat problems where the semiclassical trajectories have different energies, 
as is necessary in calculations involving time delay (this program has already been 
carried out for broken time reversal symmetry in El)- Another possibility is the 
treatment of the proximity gap in Andreev billiards 1501 . 

Similar ideas might also be applied to closed systems, allowing the calculation of 
spectral correlation functions and providing justification for the celebrated Bohigas- 
Giannoni-Schmit conjecture EUEl that they are all described by RMT. 
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Appendix 

Consider a (M + 1)—dimensional orthogonal matrix O , and let S be its A-dimensional 
upper-left square corner. It was shown in [53] that, when O is uniformly distributed in 
the orthogonal group 0(M+ 1) with Haar measure, S becomes a random matrix whose 
distribution is given by P(S) = det(l — S'S' T ) M / 2-JV , where the normalization 
constant is Z T = f dS det(l — SS T ) M / 2 ~ N . 

Suppose we wish to compute the average value of lit zl Oa k b k , with 1 <a kl h k < 
N . Since all the elements belong to S, it is obvious that 

( 2 n \ r rlQ 2n 

II ) = f- det(l - SS T ) M ' 2 - N n S akbk . (45) 

fc=1 / 0(M+ 1) T k =1 

This is an exact relation between the statistics properties of orthogonal matrices and 
those of general real matrices. 

Introduce now the different normalization constant Z = J dS'e _ ~ Tr,ss , which 
is the one required in the use of Wick’s rule. As discussed in the text, we wish to 
compute the limit as N —> 0 of Eq. (l45l) , in the form 

2 77, 

Jr™ a ^ d6t(1 " SS T ) M ' 2 - N J] S akbk • (46) 

J k=1 

It is not difficult to show that limjv-^o -yy = 1 ; and we therefore arrive at Eg. (1551) . 
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